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We study the probability distribution of the pseudo-critical temperature in a mean-field and in 
a short-range spin-glass model: the Sherrington-Kirkpatrick (SK) and the Edwards- Anderson (EA) 
model. In both cases, we put in evidence the underlying connection between the fluctuations of 
the pseudo-critical point and and the Extreme Value Statistics of random variables. For the SK 
model, both with Gaussian and binary couplings, the distribution of the pseudo-critical temperature 
is found to be the Tracy- Widom distribution. For the EA model, the distribution is found to be 
the Gumbel distribution. Being the EA model representative of uniaxial magnetic materials with 
quenched disorder like Feo.sMno.sTiOa or Euo.sBao.sMnOs, its pseudo-critical point distribution 
should be a priori experimentally accessible. 
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Disordered uniaxial magnetic materials having a glassy 
behavior like Fe .5Mn .5TiO 3 [1] and Euo. 5 Bao.5Mn0 3 [2] 
have interested physicists for decades. Since the first pi- 
oneering work of Edwards and Anderson (EA) [3] , these 
systems have been studied by means of spin-glass models 
with quenched disorder, which were later considered in 
their mean-field version by Sherrington and Kirkpatrick 
(SK) [4]. In the thermodynamic limit, Parisi's solution 
for the SK model [5] predicts a phase transition at a finite 
critical temperature separating a high-temperature para- 
magnetic phase from a low-temperature glassy phase. 
Differently, for the EA model there is no analytical so- 
lution and the existence of a finite-temperature phase 
transition relies entirely on numerical simulations [6]. 

Even though criticality in a physical system can emerge 
only in the thermodynamic limit [7, 8], in laboratory and 
numerical experiments the system size is always finite: 
singularities of physical observables are smeared out and 
replaced by smooth maxima. In order to characterize 
the critical point of finite-size systems, a suitably-defined 
pseudo-critical temperature must be introduced, e.g. the 
temperature at which such maxima occur. In finite-size 
systems with quenched disorder such a pseudo-critical 
temperature is a random variable depending on the real- 
ization of the disorder. The characterization of the dis- 
tribution of the pseudo-critical point and of its scaling 
properties is still an open problem which draw the atten- 
tion of physicists since the very first works of Harris [9- 
13]. Further studies of such distributions in spin glasses 
have been performed in a recent work [14], where some 
of the authors showed a connection between the fluctua- 
tions of the pseudo-critical temperature of the SK model 
and the theory of of Extreme Value Statistics (EVS) of 
correlated random variables. 

The EVS of Independent Identically Distributed (IID) 
random variables is a well-established problem: a fun- 
damental result [15] states that the limiting Probabil- 
ity Distribution Function (PDF) of the maximum of 
IID random variables belongs to three families of dis- 



tributions: the Gumbel, Frechet or Weibull distribution. 
Much less is known about EVS of correlated random vari- 
ables. A noteworthy case of an EVS distribution of corre- 
lated random variables that has been recently discovered 
is the Tracy- Widom (TW) distribution [16], describing 
the fluctuations of the largest eigenvalue of a Gaussian 
random matrix. The TW distribution has been found 
to describe the fluctuations of observables of a broad 
number of physical and mathematical models, like the 
longest common sequence in a random permutation [17], 
directed polymers in disordered media [18] and polynu- 
clear growth models [19], which can be described by the 
Kardar-Parisi-Zhang equation [20, 21]. Recently the TW 
distribution has been found to describe the conductance 
fluctuations in two- and three- dimensional Anderson in- 
sulators [22, 23] and measured in growing interfaces of 
liquid-crystal turbulence [24, 25] experiments. 

In this Letter we study the distribution of the pseudo- 
critical temperature in the SK and in the EA model by 
means of numerical simulations. Our numerical findings 
show that the fluctuations of the pseudo-cirtical temper- 
ature of the SK model both with Gaussian and binary 
couplings are described by the TW distribution. This 
result suggests that the features of the fluctuations of 
the pseudo-critical temperature are universal, i. e. sta- 
ble with respect to the distribution of the disorder. To 
our knowledge, this is the first time that the ubiquitous 
TW distribution is shown to play a role in spin glasses. 
Moreover, our numerical analysis shows that the fluctu- 
ations of the pseudo-critical point of the EA model are 
described by the Gumbel distribution. These two results 
shade light on the role played by EVS in spin glasses. 

To pose the problem, let us consider a system of N 
spins Si = ±1 located at the vertices of a graph, in- 
teracting via the Hamiltonian H[S] = —J2u j) JijSiSj, 
where the sum runs over the interacting spin pairs (i,j). 
For the SK model with Gaussian couplings (GSK) and 
for the SK model with binary couplings (BSK) the inter- 
acting spin pairs are all the distinct pairs. The couplings 



Jij are IID Gaussian random variables with zero mean 
and variance 1/N for the GSK model [4], and are equal 
to ±l/\/~N with equal probability for the BSK model 
[26] . For the EA model the interacting spin pairs are the 
nearest-neighbors pairs on a three-dimensional cubic lat- 
tice with periodic boundary conditions, and Jy are IID 
random variables equal to ±1 with equal probability [3]. 
For the BSK and EA model, the binary structure of the 
couplings allowed for the use of an efficient asynchronous 
multispin-coding simulation technique [6] , yielding an ex- 
tensive number of disorder samples and system sizes. 

Let us now define the physical observables used to 
carry on the numerical analysis of the problem. Given 
two real spins replicas S l ,S 2 , their mutual overlap q = 
jjJ2i=i^i^f i s a physical quantity characterizing the 
spin-glass transition in the thermodynamic limit [5, 6]: 
(q 2 )Afi) = if fi < fic, (q 2 )j(fi) > if fi > fi c , where 
(• • • ) j denotes the thermal average performed with the 

Boltzmann weight defined by the Hamiltonian H[S], 
fi = 1/T is the inverse temperature, and 7TT stands for 
the average over quenched disorder J = {Jij}ij. The 
finite-size inverse pseudo-critical temperature fi c j of a 
sample with a realization J of the disorder can be de- 
fined as the value of fi at which {q 2 )j{fi) significantly 
differs from zero, i. e. becomes critical. This qualitative 
definition is made quantitative by setting 



(q 2 )j(Pcj) = (q 2 )j(f3?) 



(1) 



Both for the GSK and BSK model, fi^ is chosen to 
be the average critical temperature at size N, which is 
defined as the temperature at which the Binder ratio 

2 

B = 1/2(3 — (q A )j I '(q 2 )j ) of a system of size N equals 
the Binder ratio of a system of size 2N. For the EA model 
we simply take fi^ to be equal to the infinite-size criti- 
cal temperature f3 c = 0.855 [27], because in this case the 
Binder ratios cross at a temperature which is very close 
to the infinite-size critical temperature (3 C . The definition 
(1) and fi^ are qualitatively depicted in Fig. 1. The dis- 
tribution of fie j can be characterized by its mean fi c j, 

2 

its variance a\ N = fi 2 j — fi c j and by the PDF pn{xj) 



cj' 

of the natural scaling variable xj = (fi c j — ficj)/o~pN- 
We can expect that, to leading order in N, apN ~ N^^ 
and that for large N pn(xj) converges to a nontrivial 
limiting PDF Poo(xj). 

Sherrington- Kirkpatrick Model - Let us start dis- 
cussing the distribution of fi c j for the GSK and 
BSK model. Monte Carlo (MC) simulations have 
been performed with Parallel Tempering (PT) for sys- 
tem sizes N = 32,64,128,256 (GSK) and N = 
16,32,64,128,256,512,1024,2048,4096 (BSK), allowing 
for a numerical computation of (q 2 )j and so of fi c j 
for several samples J . The data show that as the sys- 
tem size N is increased, fi c j approaches fi c . Setting 
T cJ = l/fi cJl a 2 TN =Tf^~T^ 2 ~ N-*, the power 
law fit of (Xt n shown in Fig. 2 gives the value of the scal- 
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FIG. 1: Square value of the overlap (q 2 )j for a sample J 
(dashed curve) for the BSK model with N = 128, its average 
{q 2 )j over the samples J (solid curve) as a function of the 
inverse-temperature /3, and critical temperatures P^ and fic j. 
The dashed vertical lines depict the definition (1) of fi c j- 
Inset: Binder parameter B as a function of the temperature 
T for N = 128, 256 and average pseudo-critical temperature 
T c 128 , with T? = l//3f. 



ing exponent = 0.31 ± 0.07 (GSK) and (j) = 0.34 ± 0.05 
(BSK). These values of <fi are both consistent with the 
value = 1/3 one would expect from scaling arguments 
by considering the variable y = iV 1 / 3 ^ - T c ) [28]. 

The PDF pn of the rescaled variable xj is depicted 
in Fig. 2. The curves pn(xj) collapse quite satisfy ingly 
indicating that we are close to the asymptotic regime 
N — > oo. Even though one could naively expect the 
fluctuations of the pseudo-critical point to be Gaussian, 
Fig. 2 shows that this is not the case. 

To understand this fact, let us recall the analysis pro- 
posed in a recent work [14] by some of the authors. 
In order to study the sample-to-sample fluctuations of 
the pseudo-critical temperature one uses the Thouless- 
Anderson-Palmer approach for the SK model. In the 
TAP approach a free energy function of the local mag- 
netization is built up for any sample J of the disorder, 
and its Hessian matrix Hij calculated at the paramag- 
netic minimum is a random matrix in the GOE ensem- 
ble. In the thermodynamic limit, the spectrum of H^ is 
described by the Wigner semicircle, centered in 1 + fi 2 
and with radius 2/3. The critical temperature fi c = 1 
of the SK model is identified as the value of fi such 
that the minimal eigenvalue of Hij vanishes. In [14] 
the fluctuations of the pseudo-critical temperature are 
investigated in terms of the fluctuations of the minimal 
eigenvalue of H^ . One introduces a definition of pseudo- 
critical temperature fi c j, which is different from that 
considered in the present work. The finite-size fluctua- 
tions of ficj are found to be described by the relation 
fi c j = fic — xj-/(2iV 2 / 3 ), where \J i s distributed ac- 
cording to the TW distribution in the high-temperature 
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FIG. 2: Distribution of the pseudo-critical point both for 
the SK model with Gaussian couplings (GSK) and for the 
SK model with binary couplings (BSK). PDF pn(xj) of 
the rescaled critical temperature xj for systems sizes N — 
32,64,128,256 with 1.6 x 10 4 < S < 4.7 x 10 4 (GSK) and 
2.9 x 10 4 < S < 9.8 x 10 4 (BSK) disorder samples, Tracy- 
Widom distribution Ptw(xj) (solid curve) and Gaussian dis- 
tribution pa{xj) (dashed curve), both with zero mean and 
unit variance. The plot has no adjustable parameters, and is 
in logarithmic scale to highlight the behavior of the distribu- 
tions on the tails. Top inset: width cttjv for the BSK as a 
function of JV and fitting function /(JV) = aN^ 4, + bN~ 2 ' 1 ' , 
yielding <fi — 0.34 ± 0.05. Bottom inset: same plot as in the 
main plot in linear scale. 



region j3 c j < 1. According to Fig. 2, MC simulations 
confirm this analysis: the limiting distribution of pn {x j) 
is described with good accuracy by the TW distribution 
in the high-temperature regime j3 c j < 1 (xj < 0). The 
TW distribution is robust with respect to the choice of 
the disorder distribution and to the definition of pseudo- 
critical temperature. On the other hand, since the expo- 
nent </> obtained from MC simulations is not compatible 
the exponent 2/3 of /3 c j, we conclude that the scaling 
exponent is definition-dependent [14, 29]. 

Edwards-Anderson Model - The same analysis has been 
performed for the three-dimensional EA model. Physi- 
cal observables have been computed with PT for system 
sizes N = L 3 with L = 4, 8, 12, 16. Similarly to the SK 
model, the width ap n of the distribution of the pseudo- 
critical point (3 C j shrinks to zero as the system size N 
is increased: a power law fit apN — aN~^ gives the 
value of the scaling exponent <fr = 0.23 ± 0.03 (inset of 
Fig. 3). The PDFs pn{xj) of the rescaled critical tem- 
perature seem to have a finite limit as N is increased, as 
depicted in Fig. 3, and this limit coincides with the Gum- 
bel distribution. Both (j> and the PDF have the following 
interesting features. As far as the exponent </> is con- 
cerned, we recall [10] that for systems known to be gov- 





/4> 

/.« 

7 1^ 
' fi 

fo "5 
f ■& 

; o 




' 








1 




PTWfe) ■ 
PGufe) ■ 


0.25 
0.2 

gi 0.15 


-X x "an X - ' 

x x 

'"x --. 
- /w - x : 


0.1 










□ 
o 

V 

♦ 


0.01 


1 
N 
N 
0.6 

0.4 

0.2 




= 12 3 

= 16 3 


% 


64 512 4096- 

N 






»\. 


- 


\ \ '•. O 


0.001 


-2 -1 


12 3 
X J 


4 5 




\ \ w 



-2-1012345 

FIG. 3: Distribution of the pseudo-critical point for the 
EA model. PDF pn(xj) of the rescaled critical temperature 
xj for systems sizes N = 4 3 , 8 3 , 12 3 , 16 3 with 2.4 x 10 4 < 
S < 3.2 x 10 4 disorder samples, Gaussian distribution pg(xj) 
(solid curve), Tracy- Widom distribution Ptw(xj) (dashed 
curve) and Gumbel distribution PGu{xj) (dotted curve), all 
with zero mean and unit variance. The plot has no adjustable 
parameters, and is in logarithmic scale to highlight the behav- 
ior of the distributions on the tails. Top inset: width opN as 
a function of JV, and fitting function /(JV) = aN~^, with 
scaling exponent <j> — 0.23 ± 0.03 . Bottom inset: same plot 
as in the main plot in linear scale. 



erned by a random fixed point like the EA model it was 
predicted that the scaling exponent satisfies l/<p = dv, 
where d is the dimensionality of the system. The value 
of the critical exponent v = 1.8 ±0.2 for the EA model is 
known from numerical simulations [6] , yielding a value of 
4> = 0.19 ± 0.02 which is compatible with that measured 
from the fluctuations of the critical temperature. As far 
as the limiting distribution Poo^j) is concerned, we re- 
call that [30] a disordered system like the EA behaves 
as an ensemble of independent sub-systems S\ , . . . , Sm , 
where each sub-system Si has a random local critical tem- 
perature /3£, the local critical temperatures {/3*}i being 
IID random variables depending on the local realization 
of the disorder. We can argue that, for a single realiza- 
tion of the disorder J ^ the pseudo-critical temperature 
j3 c j results from the fact that j3 has to be taken large 
enough to bring all of the sub-systems {Si}i to critical- 
ity. Thus, /3 C j is the maximum over the ensemble of the 
local critical temperatures fi c j = max^ (3 l c . If this picture 
is correct, j3 c j is distributed according to one of the the 
EVS limiting distributions of independent variables [15]: 
the Gumbel, Frechet, or Weibull distribution. Assuming 
that the distribution of /3* decays exponentially for large 
(3 C , the distribution of /3 c j is the Gumbel one. We want 
to stress that this argument would not hold for the SK 
model, where there is no geometric structure. 

Conclusions - In this Letter we have performed a 
numerical analysis of the distribution of the pseudo- 



crtical temperature in two mean-field spin glasses, the 
Sherrington-Kirkpatrick model with Gaussian couplings 
(GSK) and with binary couplings (BSK), and in a short- 
range spin glass, the Edwards-Anderson (EA) model. 
The analysis for the GSK and BSK model shows that 
the distribution of the pseudo-critical temperature in the 
high temperature phase is described with good accuracy 
by the Tracy- Widom (TW) distribution, as suggested by 
an analytical prediction previously published by some of 
the authors [14]. To our knowledge, this is the first time 
that the TW distribution is shown to play a role in spin 
glasses. The fact that both the GSK and BSK yield the 
TW distribution suggests that the TW distribution is 
universal with respect to the bonds' distribution. 

The analysis pursued for the three-dimensional EA 
model shows that the liming distribution of the pseudo- 
critical temperature is the Gumbel distribution. An ar- 
gument to understand this result has been proposed. 
These two numerical analyses put in evidence a connec- 
tion between the critical regime of spin-glass models and 
the Extreme Value Statistics theory which has never been 
proposed heretofore. 

The present Letter opens several perspectives. As far 
as the SK model is concerned, we recall that the TW 
distribution describes typical fluctuations of the maximal 
eigenvalue of a Gaussian Orthogonal Ensemble random 
matrix, while the large deviations regime of these fluc- 
tuations has been studied only recently [31]. It would 
be interesting to study numerically the large deviations 
regime of the fluctuations of the critical temperature, 
where the distribution of the pseudo-critical point could 
be described by the large deviations function derived in 
[31]. It would be also interesting to consider the case 
where the couplings J,j are Gaussian with a positive bias 
Jo [32]. Depending on the value of Jq, the SK model has 
a phase transition from a paramagnetic to a spin-glass 
phase or from a ferromagnetic to a mixed phase [32]: 
it would be interesting to investigate both analytically 
and numerically the fluctuations of these pseudo-critical 
points. Moreover, in order to bridge the gap between 
a mean-field and a short-range interactions regime, it 
could be interesting to investigate the fluctuations of the 
pseudo-critical temeperature in spin-glass models with 
tunable long-range interactions, like that introduced in 
[33]. As far as the EA model is concerned, it would 
be interesting to test experimentally the scenario found 
here in Feo.sMno.sTiOs [1] or Euo.sBao.sMnOa [2] spin- 
glass materials. Indeed, ac-susceptibility measurements 
in these systems show [1] that the spin-glass critical tem- 
perature can be identified as the temperature where the 
susceptibility has a cusp. Accordingly, the pseudo-critical 
point could be identified and measured, and one could 
test whether the resulting rescaled pseudo-critical point 
distribution converges to the Gumbel distribution as the 
system size is increased. 
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